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A piezoelectric mechanism of orientation of stripes in two-dimensional quantum Hall systems in 
GaAs heterostructures is considered. The anisotropy of the elastic moduli and the boundary of the 
sample are taken into account. It is found that in the average the stripes line up with the [110] axis. 
In double layer systems the wave vector of the stripe structure rotates from the [110] to [100] axis if 
the period of density modulation becomes large than the interlayer distance. From the experimental 

. point of view it means that in double layer systems anisotropic part of resistivity changes its sign 

' under variation of the external magnetic field. 

I. INTRODUCTION 

A homogeneous state of two-dimensional electron gas is unstable at low densities and low temperatures. Under such 

• conditions a transition into a Wigner crystal state occurs. The minimum of the energy corresponds to the triangular 
lattice [jl] for the case of classical Wigner crystal. Recently, non-uniform electron states in quantum Hall systems have 
become a subject of great attention. One can expect that many new states with a spatial modulation of the electron 
density are realized in these systems. In quantum Hall ferromagnets(QHF) skyrmions may form lattice structures 
H (skyrmions in QHF obey the electric charge). Skyrmions are spatially extended objects and the square lattice is 
preferable at moderate densities. Among recent experimental results the most intriguing one is the observation of 
strong anisotropy of resistivity at filling factors v = N + 1/2 (N is the integer number, N > 4) ||Q. The effect is a 

£Zj ' signature of a stripe phase state of electrons of the half-filled Landay level [g|fj] . 

It has been found in experiments |3|,^| that the low-resistivity direction is alinged with the [110] axis. Therefore there 
is an anisotropic electron-electron interaction causes preferable orientation of the stripes relative to the crystallography 
axes of the host material. Such an interaction should be comparatively weak, because the effect of rotation of the 
stripes in tilted magnetic fields has been observed @^]. 

The modulated electron structures emerge as a result of competition between the direct Coulomb and exchange 
interaction. In systems with cubic symmetry these interactions are isotropic. They cannot cause the stripe to line 
up with a particular axis. GaAs is the piezoelectric crystal. Therefore one can suppose that piezoelectric interaction 
\ (which remains anisotropic in cubic crystals) plays the main role in orientation of stripe structures. 

The anisotropy of the electron-electron interaction in piezoelectrics was studied previously by Rashba and Sherman 
in U . They investigate the influence of piezoelectric interaction on symmetry of Wigner crystal states. The isotropic 
model for the elastic media was used in Ref. ||. The results B cannot be applied directly to GaAs where the 
anisotropy of the elastic moduli is rather strong. 

In this paper we consider the piezoelectric mechanism of orientation of the stripe structures in GaAs. We take into 
account the anisotropy of the elastic moduli and the influence of the surface of the sample. The major conclusion is 
that if the electron layer is parallel to the (001) plane the energy of the stripe phase is minimal when the angle <f) 
between the wave vector of the stripe structure and the [100] axis is in the interval 30° -j- 60°. The potential curve as 
the function of <j) shows a rather shallow shape in this interval. The average direction of the wave vector corresponds 
O ■ to the orientation observed experimentally. 

Double layer stripe structures are also considered in the paper. We find that the stripes are aligned with the [100] 
• i-H , or [010] axis, if the period of the density modulation becomes greater than the interlayer distance. The effect can 
be easily checked experimentally, because the period of the stripe structure, which is proportional to the magnetic 
length, is increased under lowering of the external magnetic field. 
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II. MONOLAYER SYSTEM 



At the beginning we consider the electron layer in an infinite piezoelectric cubic crystal. The properties of the host 
system are described with three elastic moduli en, C12, C44, the dielectric constant e and one piezoelectric modulus 
ei4. Here and below only the case of a layer laying in the (001) plane is analysed. The stripe structure is described 
as a charge density wave(CDW) with the wave vector b directed at the angle 4> to the [100] axis: 
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p(r) = p sm(hr p i)S(z) (1) 

where r p i is the projection of r on the (001) plane. The z axis is along the [001] direction. The density of the energy 
of the system is 

ED (T tk u ik 



where 



is the dielectric displacement, 



A = e-Bi - knpiuuu (3) 



— XiklmUlm + PlikEl (4) 



the stress tensor, Uik, the strain tensor, E, the electric field, Xikim the elastic moduli tensor, f3iki, the piezoelectric 
moduli tensor (in cubic crystals Piki = eu/2 at i =/= k ^= I and fliki — otherwise). The quantities D and er^ satisfy 
the equations 

divD = Airp 

^ = (5) 

Taking into account Eqs.(^) and boundary condition for <7j& (cr^nj, = at the free surface, n is the vector perpendicular 
to the surface) we obtain for the total energy 

E = J d 3 rF = ~ J d 3 rp(r)(p(r) (6) 

where <p is the electric potential (E = — V</?). The value of tp can be found from Eqs. (||). Rewriting Eqs.(|^) in terms 
of the Fourier-components of tp and the displacement u we obtain 

M lk V k = Qi, (7) 

where 

«=(_| + f) (8) 



u iq i=l,2,3 ^ _ f i=l,2,3 



with Ajfc — Xiklrnqilm, Ij — —^iklQkQh 

(p q , the Fourier component of the electron density). The solution of Eq.(0) is 

= M 44 (%/ > fe)/°q ( 10 ) 

(q p ; is the projection of q on the (001) plane)). Inserting tp into Eq.(@) yields 



E = ^J dq z M- 4 \h,q z ) (11) 

(S is the area of the layer). Taking into account the smallness of the piezoelectric constant, the energy can be written 
as the sum 

E = E C + E° pe + E™ (12) 



where 



c ~ ~2dT (13) 
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is the Coulomb energy, E®™, the isotropic and anisotropic parts of the energy of piezoelectric interaction between 
electrons. The anisotropic part has the form 

EH = X E C F{4>) (14) 

where \ — ef^/ecu is the small parameter of the expansion. The function F has the amplitude of order of unity. It 
depends on the ratio between the elastic moduli. 

If C12 = en — 2C44 (the case of the isotropic elastic medium) 



F(4>) = A cos 40 (15) 



where 



The answer Eq. ( |16|) was obtained in [g| with another method. Using the values of en, C44 for GaAs in Eq. ([L6|) we 
obtain A w 0.3. The minimum of the energy reaches at <j> = tt /4. But if one can use the average values of the square 
velocities of the longitudinal and the transverse sound instead of cu, C44 correspondingly, the value of A approaches 
to zero. It means that even the sign of the anisotropy remains unsettled. Thus it is a principal question to take into 
account the anisotropy of the elastic moduli. 

The integral in Eq.(|Tl|) has been computed numerically for the anisotropic case. The values of the elastic moduli 
for GaAs (en = 12.3, C12 = 5.7, C44 = 6.0, all in 10 11 dyne/cm 2 ) are used. The result for F(<p) is shown in Fig. |l| 
(curve 1). We see from the dependence obtained that the minimum of the potential energy corresponds to (j) « 30°. 
But the energy variation in the interval 30° < <f> < 60° is very small and all the configuration in this interval are 
realized with an almost equal probability. In the average the wave vector of the stripe structure is directed along the 
[110] axis. The absolute value of the anisotropic part of the energy is determined by the parameter %■ F° r GaAs 
(ei4 = 0.15 C/m 2 , e = 12.5) x ~ 2 • 10~ 4 . 




FIG. 1. The dependence of the energy (in units of x^c) on <j> (in degrees) (0 is the angle between the wave vector of 
the stripe structure and the [100] direction). Curve 1 corresponds to the infinite system, curves 2,3,4 - the systems with the 
boundary. 2- d/a = 0; 3 - d/a — 0.15; 4 - d/a = 0.5 (d is the distance between the layer and the surface of the sample, a, the 
period of the stripe structure). 
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An important question is how the result is changed under small variation of the elastic moduli. We compute the 
function F(cj>) versus the C12 modulus. If the value of c\i is decreased the local maximum at = 7r/4 transforms to 
the global minimum (at en ~ 5 • 10 11 dyne/cm 2 ). In the opposite case (increasing of this modulus) the minimum at 
(j> w 30° becomes dipper. The situation in GaAs is just in between these two cases. 

The dependence of the energy of interaction on the angle <f> remains unchanged for the stripe structure with density 
modulation of general form. If one replaces Eq. ([[]) with the sum over the wave vectors q„ = nb (n is integer) , the 
quantities E c and are multiplied with the same factor and the function F((f>) is not modified. 

In case of the electron crystal with square lattice the energy can be written as a sum over the reciprocal lattice 
vectors and each term in the sum depends on the direction of the corresponding wave vector. Therefore a concrete 
form of the density distribution should be specified to compute the energy. If the electron crystal is modelled as two 
CDW with perpendicular wave vectors the preferable orientation and the anisotropic energy are obviously the same 
as for the stripe structure. The anisotropic energy is much smaller for the triangular lattice (modelled as the sum 
of three CDW). The absolute value of the anisotropic energy is in w 10 2 times smaller than in the stripe state. The 
minimum reaches when the angle between the [100] axis and one of the CDW wave vectors is kir/6 (k is the integer 
number). Note that anisotropic contribution into the total energy of the triangular lattice does not emerge in the 
isotropic elastic medium (see Eq.(fl5|)). 

A finite thickness of the electron layer can be taken into consideration in the approach used. The corresponding 
form factor should be added into Eq. (|Ti"|). However, since the thickness of the layer is much smaller than the period 
of the stripe structure, it does not result in essential changes. 

Usually the electron layers in heterostructures are situated close to the surface of the sample ( a typical distance 
d between the layer and the surface is w 5 ■ 10 3 A). Therefore a principal question is how the influence of the surface 
may change the result obtained. To answer the question the system (||) should be solved with taking into account the 
boundary conditions. 

The solution of (|5|) can be written in the form 

Ui = Ui(z)e lhrpl + c.c. 
Lp = (^(z)e ibr ^ + c.c. (17) 

where Ui(z) and f(z) satisfy the system of differential equations 

(c 4 4(<9 2 - by) - cnb 2 x )u x - (C12 + C44)b x {b v u v - id z u z ) - ie u byd z (p = 
(c 4 4(<9 2 - b 2 x ) - cnby)u y - (C12 + c iA )by(b x u x - \d z u z ) - ieub x d z ip = 
(cn<9 2 - c 4 4& 2 )wz + i(ci 2 + c^)d z {b x u x + byUy) + e 14 ,b x b y ip = 

e(d z - b 2 )ip + 47rei 4 (i5 z (6 K u a + b y u x ) - b x b y u z ) = (18) 



with the boundary conditions 



InEq. @ 



&iz\z=d-0 = 
<f\z=d-0 = f\z=d+0 

D z \z=d-o — {d z tp)\z=d+o 

Viz\z = -n = fiz|z=+0 

<f\z=-a = <p\z=+o 
Ui\ z= -o = Ui\ z=+ o 
F> z \z=+o — D z \z=-o = — 27rip 



a x( y )z = CAA{d z u x (y) + \b x ( y )U z ) - i-^b y ( x )(p 



(19) 
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Ozz = cnd z u z + \c^{b x u x + byUy) 
D z = -ed z ip - i2nei 4: (b x Uy + b y u x ) (20) 

An analytical solution of ([l8]) can be written if the roots of the characteristic equation are found. The roots were 
computed numerically. The energy has been separated similar to Eq. (|l2|) , with the Coulomb part 

2eb e + 1 

and the anisotropic part 
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E a P : = xKm (22) 

The dependences F(4>) obtained are shown in Fig. [j] at several values of the parameter d/a (a — 2ir /b is the period 
of the stripe structure). At d/a — (the electron layer is situated on the surface of the sample) the dependence (curve 
2) is close to the one obtained for the infinite medium (curve 1). Under increasing of the ratio d/a the minimum at 
(\> w 30° becomes dipper (curve 3). But the shape of the potential remains shallow in the interval 30° < <fi < 60°. Then 
the reverse changes begins and at d/a « 0.5 the two wall structure of the potential near <p = 7r /4 disappears (curve 
4). The anisotropy is maximal in the last case. The dependence approaches the curve 1 under further increasing of 
the parameter d/a Thus, the surface of the sample does not influence essentially on the orientation of the stripes. The 
results obtained show that the piezoelectric mechanism describes the experimental situation correctly. 



III. DOUBLE LAYER SYSTEM 



In this section we consider the piezoelectric mechanism of orientation of the stripe structures in double layer systems. 
Double layer electron systems are widely used in experiments. Therefore it is of interest to extend our consideration 
to the case of two identical stripe structures in the adjacent layers. There is another more important reason to address 
this problem. In the double layer system there is an additional parameter - the ratio between the interlayer distance 
and the period of the stripe structure. The period is proportional to the magnetic length. Therefore this parameter 
can be easily changed experimentally. One can hope that the orientation of the stripes depends on this parameter. If 
it is the case, the dependence can be extracted from the resistivity data. Comparing the experimental data with the 
theoretical prediction one can verify the model under consideration. 

Basing on the results of the previous section we consider only the infinite system. In the double layer system the 
Coulomb interaction results in a relative shift of CDW in the adjacent layers. The electron density reads as 

P (r) = po sm(br pl )(S(z - |) - 5(z + |)) (23) 

(s is the interlayer distance). Following the steps of the previous sections we find 

o 2 S f 00 

E = j dq z M^(b, q z )(l - cos(q z s)) (24) 
If the elastic moduli are isotropic 

= X E C A cos 40 (25) 

where 

A = 2 r (fe i^M ( £n_ A8 £n + 

J-co (1 + z ) c 44 C 44 

Calculating the integral in Eq.(p6f) yields 

A = %{1 - ^ - e-»((l + sb)(l - £1) + (sbf^ - ^(1 - ££)))) (27) 
Id 3C44 3C44 3C44 3 C44 

The dependence of A on the parameter s/a at C11/C44 = 12.3/6 is shown in Fig. || One can see that at s/a < 1 the 
anisotropic part of the energy changes the sign and the stripes line up with the [100] or [010] direction. 

Anisotropic systems show similar behavior. For GaAs the dependence of the energy on the angle <j> is shown in Fig. 
|§| at several values of s/a. The angle cj> which corresponds to the state with minimal energy versus the parameter s/a 
is shown in Fig. |]. The energy gain of the lowest energy state comparing to the (j) — (curve 1) and <p — 7r /4 states 
(curve 2) is also presented in Fig. [|. One can see from Figs. [| || the interaction between the layers does not change 
the orientation of the stripes at s/a > 1.5. At 0.8 < s/a < 1.5 the stripe phases with 4> ^ 0, 7r/4 can be detected. At 
s/a < 0.8 the stripes line up with the [100] or [010] direction. 

The effect obtained is of great importance. The parameter a depends on the magnetic length and the stripe phases 
at v = N + 1/2 have different periods at different N. In the double layer systems the stripes will have different 
orientation versus the filling factor (along [110] at smaller N and along [100] at larger N). The prediction can be 
verified experimentally. The effect if observed will be the experimental evidence of the major role of the piezoelectric 
mechanism of orientation of the stripe structures. It can also be used for indirect measurements of the periods of the 
stripe structures. 
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FIG. 2. The dependence of the amplitude of the anisotropic energy on the interlayer distance for the double layer system in 
an isotropic elastic medium 




FIG. 3. The dependence of the anisotropic energy on the stripe direction for the double layer system in GaAs (solid curve, 
s/a=0.75; dashed curve, s/a = 1; dotted curve, s/a — 2) 
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FIG. 4. Orientation and the anisotropic energy versus the interlayer distance in GaAs. Solid curve, <f> (in degrees); dashed 
curve, the energy gain (in relative units) relative the <f> — phase; dotted curve, the energy gain relative the <fr = tt/4 phase. 



IV. CONCLUSIONS 

We have shown that the piezoelectric interaction may play an essential role in orientation of non-uniform two- 
dimensional electron structures relative to the crystallography axes of the host material. 

The anisotropic part of the energy of the stripe structure in the quantum Hall system in GaAs heterostructures 
has been found. If the electron layer is parallel to the (001) plane the minimum of the energy corresponds to the wave 
vector of the stripe structure directed at the angle 4> « 30° to the [100] axis. The potential energy shows a rather 
shallow shape at 30° < <fi < 60° and all the direction in this interval can be realized with almost the same probability. 
In the average the stripes line up with the [110] axis. 

The influence of the surface of the sample on the orientation of the stripes has been analysed. If the electron layer 
is situated on the surface of the sample (the distance d between the layer and the surface is much smaller than the 
period a of the stripe structure) the anisotropic part of the energy is almost the same as for the infinite system. If the 
ratio d/a is of order of unity some changes may take place. At d/a « 0.5 the local maximum at <j> = tt/4 transforms 
into the global minimum. At d/a « 0.15 the minimum at (j> « 30° becomes dipper and the orientation of stripes along 
a low symmetry direction can be observed at very low temperatures. 

The results obtained describes the experimental situation correctly. The absolute value of the anisotropic energy is 
small ( in w 10 4 times smaller than the Coulomb energy). Therefore external fields (for instance, the tangent magnetic 
field) may change the orientation (this effect observed experimentally). 

The influence of the piezoelectric interaction on the orientation of the stripes is investigated with reference to the 
double layer systems. For such systems the model predicts that at small interlayer distances the stripes change their 
direction and line up with the [100] axis. The re-orientation takes place if the interlayer distance becomes smaller 
than the period of the stripe structure. The prediction can be checked experimentally since the period is proportional 
to the magnetic length, e.g. it will change under magnetic field variation. 

If the effect is observed it can be considered as the experimental evidence of the major role of the piezoelectric 
mechanism of orientation of stripe structures. 

This work was supported in part by the INTAS grant No. 97-0972. 
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